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Abstract. A survey of currently known planet-hosting stars indicates that approximately 
25% of extrasolar planetary systems are within dual-star environments. Several of these 
systems contain stellar companions on moderately close orbits, implying that studies of 
the formation and dynamical evolution of giant and terrestrial planets, in and around bi- 
nary star systems have now found realistic grounds. With the recent launch of the space 
telescope CoRoT, and the launch of NASA's Kepler satellite in 2009, the number of such 
dynamically complex systems will soon increase and many more of their diverse and inter- 
esting dynamical characteristics will soon be discovered. It is therefore, both timely and 
necessary, to obtain a deep understanding of the history and current status of research 
on planets in binary star systems. This chapter will serve this purpose by reviewing the 
models of the formation of giant and terrestrial planets in dual-star environments, and by 
presenting results of the studies of their dynamical evolution and habitability, as well as 
the mechanisms of delivery of water and other volatiles to their terrestrial-class objects. 
In this chapter, the reader is presented with a comprehensive, yet relatively less techni- 
cal approach to the study of planets in and around binary stars, and with discussions on 
the differences between dynamical characteristics of these systems and planetary systems 
around single stars. 



1.1 Introduction 

The concept of a "world with two suns" has been of interest to astronomers for 
many years. Many scientists tried to understand whether planets could form in 
binary star systems, and whether the notion of habitability, as we know it, could 
be extended to such environments. Although as a result of their respective works, 
many dynamical features of binary-planetary systemtQ have been discovered, until 
recently, the subjects of their studies were, in large part, hypothetical. There was 
no detection of a planet in and/or around a binary system, and planet detection 
techniques had not advanced enough to successfully detect planets in dual-star 
environments. 

The discovery of extrasolar planets during the past decade has, however, changed 
this trend. Although the candidate planet-hosting stars have been routinely chosen 
to be single, or within wide (>100 AU) binaries 0, the precision radial velocity 



A binary-pla netary system is a dual-star system that also hosts planetary bodies. 
As shown by iNorwood fc Haghighipourl (j2002l l. the perturbative effect of the stellar 
companion on the dynamics of a planetary system around a star becomes important 
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Fig. 1.1. Velocity re siduals to 7 Cepliei after sub tracting a second order fit to its original 
radial velocity data dCampbe ll, Wal ker fc, Yan g I, llOSd). The residuals show periodicity 
implying the possible existence of a planetary companion. 



technique has been successful in detecting planets around the primaries of three 
moderately close (<40 AU) dual-star systems. As a result, during the past few years, 
the topic of planets in binaries, once again, found its way to the mainstream research 
and has now become a real scientific issue that demands theoretical explanations. 

The first detection of a planet in a binary system was reported by Campbell, 
Walker & Yang in 1988. In an attempt to identify planetary objects outside our 
solar system, these authors measured the variations in the radial velocities of a 
number of stars, and reported the possibility of the presence of a Jovian-type body 
around the star 7 Cephei (figure g.l. lCampbell. Walker fc Yang I. flOSSh. Th is star, 
that is a Kl IV sub-giant with a mass of 1.59 solar-masses (jFuhrmann is the 



primary of a binary system with a semimajor axis of 18.5 AU an d an eccentricity 
of 0.36 (fCriffin. Carguillat. fc Ginestet I . l2002t iHatzes et all . l2003l). The secondary 



of th i s system is an M dwarf with a mass of 0.44 solar-masses (jNeuhauser et al 



2007t [Torres 1 . 120071 ) . Initial radial velocity measurements of 7 Cephei implied that 



this star may be host to a giant planet with a probable mass of 1.7 Jupiter-mass es, 



198ah . 



in an orbit with a semimajor axis of 1.94 AU (jCampbell. Walker fc Yang 

Unfortunately, the discovery of the first binary-planetary system, which could 
have also marked the detection of the first planet outside our solar system, did not 
withstand skepticism. In an article in 1992, Walker and his colleagues attributed 
their measured variations of the radial velocity of 7 Cephei to the chromospheric 



when the binary separation is smaller than 100 AU. At the present, approximately 
25% of extrasolar planetary systems detected by radial velocity technique are in binary 
systems with separations ranging from 250 to 6000 AU. 
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activities of this star, and announced t hat the possibil i ty of t he existence of a giant 



planet around 7 Cephei may be none (jWalker et al. 1 . 119921 ). It took observers an 



additional 12 years to monitor 7 Cephei and measure its radial velocity to arrive at 
the conclusion that the previously observed variations were not due to stellar activ- 
ities and were in fact representative of a planetary companio n (figure 9.2). It wa s 



the discovery of a giant planet in the binary system of GL 86 ( Queloz et al. . 20o3), 



and the (re-)announ cement of the detection of a giant planet in 7 Cephei system 
(jHatzes et al i;hat opened a new chapter in the theoretical and observational 

studies of extrasolar planetary systems. 

The fact that giant planets exist in moderately close (< 40 AU) binary sys- 
tems has confronted dynamicists with many new challenges. Questions such as, 
how are these planets formed, can binary-planetary systems host terrestrial and/or 
habitable planets, how are habitable planets formed in such dynamically complex 
environments, and how do such planets acquire the ingredients necessary for life, 
are among major topics of research in this area. This chapter is devoted to review 
these issues and present the current status of research on the formation of planet in 
dual-star systems and habitability of terrestrial bodies in and around binary stars. 

The chapter begins with a review of the dynamics of a planet in a binary star 
system. In general, prior to constructing a theory for the formation of planets, it 
proves useful to study whether the orbit of a planet around its host star would 
be stable. In a binary system, such studies are of quite importance since in these 
systems the perturbation of the stellar companion may dictate the possibility of 
the formation of planetary bodies by affecting the stability and dynamics of smaller 
objects. 

The formation of planets in binary star systems is reviewed in the third 
section. Although the study of the dynamics of planets in and around binary 
stars dates back to approximately forty year ago, the formation of planets in 
these systems is an issue that is still unresolved. In spite the observational ev- 
idence that indicate m ajority of main and pre-main seq u ence stars are formed 
in binaries or cluster s jAbt" . 1979t Duauennov &: Mayor . 1991 : Mathieu . 1994; 



Mathieu et aTl . I2OO0I : IWhite fc Ghez I . l200l[) . and in spite the detection of po 



tentiall y planet-forming environments in and around binary stars (figure 9.3 , 
also seelMathieu l.ll994ilAkeson. Koerner fc Jensen .119981: iRodriguez et al. I . Il998t 



IWhite et al. I . Il999l : ISilbert et al. I . l200d iMathieu et al. 1 . 12000[ ). planet formation 
theories are still unclear in explaining how planets may form in multi-star envi- 
ronments. The focus of section 9.3 is on discussing the formation of giant and 
terrestrial planets in moderately close binary-planetary systems, and reviewing the 
current status of planet formation theories in this area. 

The habitability of a binary system is presented in section 9.4. The notion of 
habitability is defined based on the habitability of Earth and life, as we know it. 
Such a definition requires a habitable planet to have the capability of retaining 
liquid water in its atmosphere and on its surface. The latter is determined by the 
luminosity of the central star, the size of the planet, and also the distribution of 
water in the protoplanetary disk from which terrestrial-class objects are formed. 
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Fig. 1.2. Radial velocity measurements of the primary of 7 Cephei binary system. The 
graph on the top shows the orbital solution for the planet (solid line) and the residual 
velocity measurements after subtracting the contribution due to the binary companion 
(data points). The graph on the bottom depicts the phased residual radial velo city mea- 
surem ents (data points) compared to the planet orbital solution (solid line) (Hatzes et al.l , 
l2003l ). 
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In section 9.4, a review of the current status of the models of habitable planet 
formation in and around binary systems are presented, and their connections to 
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Fig. 1.3. Image of the circumbinary disk of GG Tau l|Krist et al. I , |2005| ) taken by Ad- 
vanced Camera for Surveys (ACS) on board of Hubble Space Telescope (HST). The binary 
separation is 35 AU. The locations of the binary components are marked with crosses. 

models of terrestrial planet formation and water- delivery around single stars are 
discussed. Finally the chapter ends by discussing the future prospects of research 
in the field of planets in binaries. 

1.2 Dynamical Evolution and Stability 

In general, one can consider orbital stability synonymous with the capability of an 
object in maintaining its orbital parameters (i.e., semimajor axis, eccentricity, and 
inclination) at all times. In other words, an object is stable if small variations in 
its orbital parameters do not progress exponentially, but instead vary sinusoidally. 
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Instability occurs when the perturbative forces create drastic changes in the time 
variations of these parameters and result in either the ejection of the object from 
the system (i.e., leaving the system's gravitational field), or its collision with other 
bodies. 




The concept of stability, as explained above, although simple, has been defined 
differently by different authors. A review by Szebehely in 19 84 lists 50 d i fferen t 
definitions for the stability of an object in a multi-body system (| Szebehely 1 . 1 19841 ). 



ty ^ „ „ „ , ^ 

For instance. [Harrington I ( 19771 ) considered the orbit of an object stable if, while 
numerically integrating the object's orbit, its se r nimaj or axi s and orbital eccentricity 
would not undergo secular changes. ISzebehelv" ( 1980[ ). and lSzebehelv h McKenziel 
( 198l[ ). on the other hand, considered integrals of motion and curves of zero velocity 
to determine the stability of a planet in and around binary stars. In this chapter, 
the values of the orbital eccentricity of an object and its semimajor axis are used 
to evaluate its stability. An object is considered stable if, for the duration of the 
integration of its orbit, the value of its orbital eccentricity stays below unity, it does 
not collide with other bodies, and does not leave the gravitational field of its host 
star. 

The study of the stability of a planetary orbit in dual stars requires a detailed 
analysis of the dynamical evolution of a three-body system. Such an analysis itself 
is dependent upon the type of the planetary orbit. In general, a pla netary-class 



object may have thr e e typ es of orbit in and around a binary star system. ISzebehelv 



(1983) and Dvorak ( 1982 ) h ave div ided these orbits into three different categories. 
As indicated bv ISzebehelv ( 1980l ). a planet may be in an inner orbit, where it 
revolves around the primary star, or it may be in a satellite orbit, where it revolves 
around the secondary star. A planet may also be revolving the entire binary sys tem 
in which case its orbit is called an outer orbit. As classified bv iDvorak" 

I (|l982h . on 

the other hand, a study of the stability of resonant periodic orbits in a restricted, 
circular, three-body system indicates that a planet may have an S-type orbit, where 
it revolves around only one of the stars of the binary, or may be in a P-type orbit. 
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where it revolves the entire binary system (figure 9.4). A planet may also be in 
an L-type orbit where it librates in a stable orbit around the L4 or L5 Lagrangian 
points. 

The rest of this section is devoted to a review of the studies of the stability of 
planets in S-type and P-type orbits. Given that in a binary star system, a planet is 
subject to the gravitational attraction of two massive bodies (i.e., the stars), it would 
be important to understand how the process of the formation of planetary objects, 
and the orbital dynamics of small bodies would be affected by the orbital charac- 
teristics of the binary's stellar components. In general, except for some special cases 
for which analytical solutions may exist, such studies require numerical integrations 
of the orbits of all the bodies in the system . In th e past, prior to the invention of 
symplectic integrators ( Wisdom &: Holman I . Il991 ). which enabled dynamicists to 
extend the studies of the stability of planetary systems to several hundred million 
years, and before the development of fast computers, majority of such studies were 
either limited to those special cases, or were carried out numerically for only a small 
numbe r of binary's orbita l period. Exainples o f suc h studies can be fou nd in arti- 
cles bv lGraziani fc Black I (|l98l[ ). lBlackl |l982h . and lPendleton fc Black I (|l983h . in 
which the authors studied the orbital stability of a planet in and around a binary 
star. By numerically integrating the equations of the motion of the planet, this 
authors showed that, when the stars of the binary have equal masse s, the orbital 



stabi lity of the planet is independent of its orbital inclination (also see [Harrington 



Their integrations also indicate that, when the mass of one of the stellar 
components is comparable to the mass of Jupiter, planetary orbits with inclinations 
higher than 50° tend to become unstable. 

The invention of symplectic integrators, in particular routines that have been 
designed specif ically for the purpose o f integrating orbits of small bodies in dual- 
star system^ ()Chambers et al. have now enabled dynamicists to extend 
studies of planet formation and stability in dual-star environments to much larger 
timescales. In the following, the results of such studies are discussed in more detail. 



1.2.1 Stability of S-type orbits 

As mentioned above, instability occurs when the perturbative effects cause the semi- 
major axis and orbital eccentricity of a planet change in such a way that either the 
object leaves the gravitational field of the system, or it collides with another body. 
For a planet in an S-type orbit, the gravitational force of the secondary star is the 
source of these perturbations. That implies, a planet at a large distance from the 
secondary, i.e., in an orbit closer to its host star, may receive less perturbation from 

Symplectic integrators, as they were originally developed bv lWisdom fc Holman I (|l99ll ), 
are not suitable for numerically integrating the orbits of small bodies in the gravitational 
fields of two massive objects. These integrators have been designed to integrate the orbits 
of planet ary or smaller objects wh en they revolve around only one massive central body. 
Recently [Chambers et al. I l|2002l ) have developed a version of a symplectic integrator 
that is capable of integrating the motion of a small object in the gravitational fields of 
two stellar bodies. 
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the bi nary companion an d may be able to sustain its dynamical state for a longer 



time ([Harrington I . Il977f) . Since the perturbative effect of the stellar companion 



varies with its mass, and the eccentricity and semimajor axis of the binary (which 
together determine the closest approach of the secondary to the planet), it is pos- 
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Fig. 1.5. Graphs of the critical semim ajor axis (ae) of an S-type b inary-planetary system, 
in units of the binary semimajor axis (|Holman fc Wiegert I . [l999l ) . The graph on the left 
shows function of the binary eccentricity for an equal-mass binary. The graph on 

the right corresponds to the variations of the critical semimajor axis of a binary with an 
eccentricity of 0.5 in term of the binary's mass-ratio. T he solid and dashed l ine o n the 
left panel depict the e mpirical formulae as reported by iHolman fc Wiegert I and 
iRabl fc Dvorak I (119881 ) . respectively 



sible to estimate an upper limit for the planet's distan ce to the star beyond w hich 
the orbit of the planet wou ld be unstable. As shown bv lRabl fc Dvorak" (1988) and 



Holman fc Wiegert I([l993), the maximum value that the semimajor axis of a planet 
in an S-type orbit can attain and still maintain its orbital stabilit y is a fimction of 
the mass-ratio and orbital elements of the binary, and is given by (jRabl fc Dvorak 
1988t iHolman fc WiegerT 



19991) 



ttc/at = (0.464 ± 0.006) + (-0.380 ± 0.010)^^ + (-0.631 ± 0.034)ef, 

-l-(0.586 ± 0.061)Aieb + (0.150 ± 0.041)eg + (-0.198 ± 0.047)^eg . (1.1) 

In this equation, Oc is critical semimajor axis , = Mi /(Mi -I- M2), Of, and Cb are 
the semimajor axis and eccentricity of the binary, and Mi and M2 are the masses of 
the primary and secondary stars, respectively. The ± signs in equation (9.1) define 
a lower and an upper value for the critical semimajor axis Oc, and set a transitional 
region that consists of a mix of stable and unstable systems. Such a dynamically gray 
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area, in which the state of a system changes from stabihty to instabihty, is known to 
exist in multi-body environments, and is a characteristic of any dynamical system. 

Equation (9.1) is an empirical formula that has been obtained by numerically 
integrating the orbit of a test particle (i.e., a massless object) at different distances 
from the primary of a binary star ( Rabl fc Dvorak I . llQSSt iHolman fc Wiegert I . 
I999I). Figure 9.5 show s this in mo re detail. Simila r stud ies have been done by 



Moriwaki fc: Nakagawa ( 2004 ). and Fatuzzo et al. ( 2006h who obtained critical 
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Fig. 1.6. Graphs of the semimajor axes (left) and eccentricities (right) of the giant planet 
(black) and binary (gr een) of 7 Caphei for different values of the eccentricities of the binary 
( HaghighipouTI. |2003| ) . The mass-ratio of the binary is 0.2. 



semimajor axes slightly larger than given by equation (9.1). 

Since the mass of a Jovian-type planet is approximately three orders of mag- 
nitude smaller than the mass of a star, such a test particle approximation yields 
results that not only are applicable to the stability of giant planets, but can also 
be us ed in identifying regions where smaller bodie s, such as terrestrial-class ob 



jects ( Quintana et al. . 20021: Quint ana fc Lissauer I . l200d Quintana et al 



200<fl 



and dust part icles ( Trilling et al. . 20071) . can have long term stable orbit^. In a 
recent article. [Trilling et al. I ( 2007 ) utilized equation (9.1) and its stability criteria 
to explain the dynamics of debris disks, and the possibility of the formation and 
existence of planetesimals in and around 22 binary star systems. By detecting in- 
frared excess of dust particles, these authors confirmed the presence of stable dust 
bands, possibly resulted from collision of planetesimals, in S-type orbits in several 
wide binaries. 



T he stability of S-type systems has been studied by many authors (iBenest 



1988lll989L[l993[ll996tlWiegert fc Holman l . ll997tlPilat-Lohinger fc Dvorak 



2002 



In applying equation (9.1) to the stability of dust particles, one has to note that this 
equation does not take into account the effects of non-gravitational forces such as gas- 
drag or radiation pressure. The motion of a dust particles can be strongly altered by 
the effects of these forces. 
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Dvorak et aTl . l2003l |2004 iPilat-Lohinger et al. 1 . 12004 iMusielak et al. 1 . 12005| ). In 



a recent article, IHaghighipour I (j2006D extended such studied to the dynamical sta- 
bility of the Jupiter-like planet of the 7 Cephei planetary system. By numerically 
integrating the orbit of this object for different values of at, Cb and (the orbital 
inclina tion of the giant planet relative to the plane of the binary), IHaghighipour I 
(|2006 h has shown that the orbit of this planet is stable for the values of the binary 
eccentricity within the range 0.2 < Cb < 0.45. Figure 9.6 shows the results of such 
integrations for a coplanar system with n = 0.2 and for different values of the binary 
eccentricity. The initial value of the semimajor axis of the binary was chosen to be 
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Fig. 1.7. Graphs of the semimajor axes (top) and eccentricities (bottom) of the giant 
planet (black) and binary (red) of 7 Caphei. The initial eccentricity of the binary at 
the beginning of num erical integration and the value of its mass-ratio were equal to 0.20 
(IHaghighipour I.I2OO6I '). 



21.5 AU. Integrations also indicated that the binary-planetary system of 7 Cephei 
becomes unstable in less than a few thousand years when the initial value of the 
binary eccentricity exceeds 0.5. 

Interesting results were obtained when the 7 Cephei system was integrated for 
different values of ip. The results indicated that for the above-mentioned range of 
orbital eccentricity, the planet maintains its orbit for all values of inclination less 
than 40°. Figure 9.7 shows the semimajor axes and orbital eccentricities of the 
system for Cb = 0.2 and for ip=5°, 10°, and 20°. For orbital inclinations larger than 
40°, the system becomes unstable in a few thousand years. 



Kozai Resonance 

An interesting dynamical phenomenon that may occur in an S-type binary, and has 
also been observed in the numerical simula tions of a few of these system s, is the 



Koza i resonance IH aghighipour , 2003, 2005j jyerrier fc Evans ..2006: Takcda fc Ra,sio 
20061 iMalmberg. Davies fc Chambers I . |2od3). As demonstrated by iKozai I (|l962l ). 



in a three-body system with two massive objects and a small body (e.g., a binary- 
planetary system) , the exchange of angular momentum between the planet and the 
secondary star, can cause the orbital eccentricity of the planet to reach high values 
at large inclinations. Averaging the equations of motion of the system over mean 



1 Extrasolar Planets in Binary Star Systems 11 



anomalies, one can show that in this case, the averaged sy stem is integrable when 



the ratio of distances are large (the Hill's approximation, iKozai 1 . 119621 ) . The La- 



grange equations of motion in this case, indicate that, to the first order of planet's 
eccentricity, the longitude of the periastron of this object, Wp, librates around a fix 
value. Figure 9.8 sh ows this for the giant plan et of 7 Cephei. As shown here, Up 



librates around 90° (|Haghighipour I . l2003l l2005l ) 



In a Kozai resonance, the longitude of periastron and t he orbital eccent ricity of 



the small body (ep)are related to its orbital inclination as (jlnnanen 1 . 1X997^ 



0.4csc^i 



(1.2) 



and 



30 p 
25 - 
20 

10 - 

5 - 

ot" 

360 p 
270 - 
^ 180 - 
90 

'I 



Semimajor Axis 



Longitude of Periastron 



— — 



27.4 54.8 

Million Years 




2.74 5.48 8.22 
Million Years 



10.96 



Fig. 1.8. Graphs of the semimajor axis and eccentricity of the giant planet (black) and 
binary (red) of 7 Cephei (top) and its longitude of periastron a n d red uced Delaunay 
momentum (bottom) in a Kozai resonance ( Hag highipour I . 12001 |2005| ). As expected, 
the longitude of the periastron of the giant planet oscillates around 90° and its reduced 
Delaunay momentum is constant. 



(el). 



1 — 5 cos(2ip) 



(1.3) 
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From equation (9.3), one can show that the Kozai resonance may occur if the 
or bital inclination of the smal l body is larger than 39.23°. For instance, as shown 
by iHaghighipourl ( 20031 2005f ). in the system of 7 Cephei, Kozai resonance occurs 
at ip = 60°. For the minimum value of ip, the maximum value of the planet's orbital 
eccentricity, as given by equation (9.4), is equal to 0.764. Figure 9.8 also shows that 
Cp stays below t his limiting value a t all times 



As shown bv iKozafl (|l962,) and llnnanen (1992), in a Kozai resonance, the dis- 
turbing function of the system, averaged over the mean anomalies, is independent 
of the longitudes of ascending nodes of the small object (the planet) and the per- 
turbing body (the stellar companion). As a result, the quantity \/ a{l — e^) cosi 
(shown as the "Reduced Delaunay Momentum" in figure 9.8) becomes a constant 
of motion. Since the eccentricity and inclination of the planet vary with time, the 




(10000 Years) 

Fig. 1.9. Graphs of the eccentricity and incl ination of the giant planet of 7 Cephei in a 
Kozai Resonance l|Haghighipour 1 12003| . [ioosh . As expected, these quantities have similar 
periodicity and are 180° out of phase. 

fact that the quantity above is a constant of motion implies that the time- variations 
of these two quantities have similar periods and, at the same time, they vary in such 
a way that when ip reaches its maximum, Cp reaches its minimum and vice versa. 
Figure 9.9 shows this clearly. 



1.2.2 Stability of P-type Orbits 



Numerical simulations have also bee n carried out for the stability of P-type orbits 



in binary-planetary systems (IZiglin . 1975: Szcbchclv & McKenzie , 1981: Dvorak , 


1984'.'l986': 'Dvorak. Froeschle. & Froeschle 1. 


1989: Kubala, Black & Szebehelv: Holman & Wiegert . 


199^ 


Broucke . 2001: Pilat-Lohinger. Funk 


Dvorak 


.,2003:lMusielak et al. . 2005). 



Similar to S-type orbits, in order for a P-type planet to be stable, it has to be at a 
safe distance from the two stars so that it would be immune from their perturbative 
effects. That is, planets at large distances from the center of mass of a binary will 
have a better chance of begin stable. This distance, however, cannot be too large 
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because at very large distances, other astronomical effects, such as galactic pertur- 
bation, and perturbations due to passing stars, can render the orbit of a planet 
unstable. 

To determine the critical value of the semimaj or axis o f a P- type planet in 
a stable orbit, preliminary attempts were made bv [Dvorak (1984), who numeri- 
cally integrated the orbit of a circumbinary planet in a circular orbit around an 
eccentric binary system and showed that planets at distances 2- 3 times the sep - 
aration of the binary have stab l e orb its. Su bsequent studies bv lPvorakl (|l986 



Dvorak. Froeschle. &: Froeschle I (|l989l ). and iHolman fc Wiegert I (|l999l ) comple- 



mented Dvorak's results of 1984 and showed that the orbit of a P-type planet will 
be stable as long as the semimajor axis of the planet stays larger than the critical 
value given by (figure 9.10) 
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Fig. 1.10. Critical semimajor a. x is as a function of the binary eccentricity in a P- 
type system (|Holman fc Wiegert I. 19991). The sq uares correspond to the result of sta- 
bility simulations by iHolman fc Wieeertl (|l999[ l and the triangles represent those of 



Dvora k. Froes chle. fc Froeschle (119891 ). The solid line corresponds to equation (9.4). As in- 



dicated by Holman fc Wiegert I l|l999i l . the figure shows that at outer regions, the stability 



of the system fades away. 
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Oc/af, = (1.60 ± 0.04) + (5.10 ± 0.05)6,, + (4.12 ± 0.09)^* 

+ (-2.22 ± 0.11)e2 + (-4.27 ± 0.l7)ebfi + (-5.09 ± 0.11)/^^ 
+ (4.61±0.36)eg/x2 . 



(1.4) 



Similar to equation (9.1), equation (9.4) represents a transitional region with a lower 
boundary below which the orbit of a P-type planet will be certainly unstable, and 
an upper boundary beyond which the orbit of the planet will be stable. The mixed 
zone between these two boundaries represents a region where a planet, depending 
on its orbital parameters, and the orbital parameters and the mass-ratio of the 
binary, may or may not be stable. Recently , by applying the sta bility criteria of 



equation (9.4) to their observational results. iTrilling et al. I (20o3) have confirmed 
the presence of stable dust band, possibly resulted from collision of planetesimal, 
around close binary star systems. 

A dynamically interesting feature of the stable region aro und the stars of a bi- 



nary i s the appearance of islands of instability. As shown by iHolman &: Wiegert 



( I999I ) islands of instability may develop beyond the inner boundary of the mixed 
zone, which correspond to the locations of (n : 1) mean-motion resonances. The ap- 
pearance of these u nstable regions have been reported by several authors under var- 



ious circumstances (Henon fc Guvot , 1970l : iDvorak . 1984 : Rabl fc Dvorak . 1988 



Dvorak. Froeschle. fc Froeschle . 1989[) . Extensive numerical simulations would be 



necessary to determine whether the overlapping of these resonances would result in 
stable P-type binary-planetary orbits. 



1.3 Planet Formation in Binaries 



Despite a wealth of articles on planets in binary star systems, the process of the 
formation of these objects is still poorly understood. The current theories of planet 
formation focus only on the formation of planets in a circumstellar disk around a 
single star, and their extensions to binary environments are limited to either the 



inner planets of our solar svstem dHeooenheimer 


. 1974, 19781 Drobvshevski . 1978: 


Diakov & Reznikov . 19801 Whitmire et al. 1. 1998 


: KortenkamtJ. Wetherill & Inabah 



or binaries resembling some of ex trasolar planets, in whic h the secondary star has a 



mass in the brown dwarf regime ( Whitmire et al. , 1998h . Although attempts have 



been ma de to extend such studies to binaries with comparable-mass stellar com- 
ponents dMarzari &: Scholl I. I2OOOI: iNelson I. [ioOoHBarbieri. Marzari fc Scholl 1 . 12002 



Quintana et al. , 2002 ; Lissauer et al. , 20041 ). the extent of the applicability of the 
results of these studies has been only to hypothetical cases since, until recently, 
there had been no observational evidence on the existence of such binary-planetary 
systems. 
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In general, it is believed that planet formation proceeds through the following 
four stages (figure 9.11): 

1) coagulation of dust particles and their growth to centimeter-sized objects, 

2) collisional growth of centimeter-sized particles to kilometer-sized bodies (plan- 
etesimals), 

3) formation of Moon- to Mars-sized protoplanets (also known as planetary 
embryos) through the collision and coalescence of planetesimals, and 

4) collisional growth of planetary embryos to terrestrial-sized objects. 

The latter is a slow process that may take a few hundred million years. During the 
first few million years of this process, at larger distances from the star, planetesimals 
and planetary embryos may form planetary cores several times more massive than 
Earth, and may proceed to form giant planets. 

In a binary star system with a moderate to small separation, the secondary 
star will have significant effects on the efficiency of each of these processes. As 
shown by Boss (2006), a binary companion can alter the structure of a planet- 
forming nebula, and create regions where the densities of the gas and dust are lo- 
cally enhanced (figure 9.12). Also, as shown bv I Artvmowicz fc Lubow" (1994), and 




Fig. 1.11. The four stages of planet formation. 



Pichardo. Sparke fc Aguilar I (|2005f ). a stellar component on an eccentric orbit can 



truncate the circumprimary disk of embryos to smaller radii and remove material 
that may be used in the formation of terrestrial planets (figure 9.13). As a result, 
it used to be believed that circumstellar disks around the stars of a binary may no t 
be massive enough to form planets. However, observations by iMathieu I (|l994f ). 
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Fig. 1.12. Structure of a circumprimary disk in a double star system. The masses of the 
primary and secondary stars are 1 and 0.09 solar-masses, respectively. The secondary star 
is at 50 AU at the top of the figure, and has an eccentricity of 0.5. The figure shows an 
area of 20 AU around the primary. The structures inside the disk have appeared after 239 
years from the beginning of the simulations. The orange st ructure on the right edge of the 
graph is an artifact of numerical simulations ||Boss I . I2OO6I ). 



Akeson. Koerner fc Jensen] ( 19981 ). iRodriguez et al. ( 1998 ). and Mathieu et al 



( 200Cl[ ) have indicated that potentially planet-forming circumstellar disks can in- 



deed exist around the stars of a binary system, implying that planet formation in 
binaries may be as common as around single stars (figure 9.14). Among these cir- 
cumstaller disks, the two well-separated disks of the system L1551 retain equivalent 
of approximately 0.03 to 0.06 solar-masses of their origin al circumstellar material s 



in a region with an outer radius of '--^lO AU (figure 9.14, [Rodriguez et al. I . Il998f ) 



The masses of these disks are compara ble to the minimum solar-mass model o f the 



primordial nebula of our solar system (jWeidenschilling 1 . 119771 : iHavashi 1 . 119811 ) , im 



plying that, planet formation in dual-star systems can begin and continue in the 
same fashion as around our Sun. 

Despite the observational evidence in support of the existence of planet-forming 
environments in moderately close binary star systems, the perturbative effect of the 
binary compan ion may not always favor planet formation. For instance, as shown by 



Nelson 



, , ■ ( 2000l). gia.nt pl anet formation cannot proceed through the disk instability 
mechanism ()Boss . 2000l ) around the primary of a binary star system with separation 
of ^ 5 AU. Also, when forming pla netary embryos, a s sho wn by iHeppenheimer I 
(|l978h . IWhitmire et al. I (|l998h . and iThebault et all (|2004l ). the perturbation of 



the secondary star may increase the relative velocities of planetesimals and cause 
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Fig. 1.13. Disli Truncation in and around binary systems (lArtvmowicz fc Lubow Lll994l l. 
The top graplis show circumstellar dislt in a binary with a mass-ratio of 0.3. Note the 
disk truncation when the eccentricity of the binary is increased from to 0.3. The bottom 
graphs show similar effect in a circumbinary disk. The mass-ratio is 0.3 and the binary 
eccentricity is 0.1. The numbers inside each graph represent the time in units of the binary 
period. The axes are in units of the binary semimajor axis. 
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Fig. 1.14. Interferometric observation of the binary system L1551 (jRodriguez et alH . 
\WM)- Two compact sources are evident in the map. The separation of the binary is 45 
AU and the disk around each core extends to approximately 10 AU. 



their collisions to result in breakage and fragmentation (figure 9.15). Results of the 
studies by these authors suggest that planetesimal accr etion will be efficient only 
in binaries with large separation [50 AU as indicated by lHeppenheimer_ dlQTSl). 26 
AU as shown by lWhitmire et al~ ( 1998h . and 100 AU as reported by Maver et al. I 
( 2005f )]. Finally, in a binary star system, the stellar companion may create unstable 



regions where the building blocks of planet s will not maintain thei r orbits and, as 



a ri 



esult, planet formation will be inhibited (jWhitmire et al. I . Il998h . 



Interestingly, despite all these difficulties, numerical simulations have shown that 
it may indeed be possible to form giant and /or terres trial planets in and around a 
dual-s tar system. Recent simulations by Boss ( 2006[ l. and Maver. Boss fc Nelson I 
(|2007l ) indicate that Jupiter-like planets can form around the primary of a binary 
star system via gravitational instability in a ma rginally unstable circumprimary 
disk (figure 9.16). On the other hand, as shown bv lThebault et al~ (2004), the core 
accretion mechanism may also be able to form giant planets around the primary of 
a binary star. However, as the results of their simulations for planet formation in 
the 7 Cephei system indicate, the semimajor axis of the final gas-giant planet may 
be smaller than its observed value. 
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Fig. 1.15. Graphs of the evolution of eccentricity (top) and encounter velocities (bottom) 
for planetesimals at th e region between 0.3 and 5 AU from the primary of 7 Cephei 
ijThebault et al. Il2004l ). The planetesimals disk in the bottom simulation was initially at 
its truncated radius of 4 AU. As shown here, the perturbative effect of the secondary star 
increases the eccentricities and relative velocities of these objects. 
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Fig. 1.16. Giant planet formation via disk instability mechanism in a binary system. The 
separation of the binary is 120 AU and it was initially on a circular orbit. The mass of 
each disk is 0.1 solar-masses. The snap shot was taken 160 years after the start of the 
simulations. Figure Courtesy of L. Mayer, A. Boss and A. Nelson . 



In regard to the formation of terrestrial planets in binary systems, in a series of 
articles, Quintana and her colleagues integrated the orbits of a few hundred Moon- to 
Mars-siz ed objects and showed that terrestrial-c l ass objects can form in and around 



binaries ( Quintana et al. , 2002t Lissauer et al. . 2004 : Quintana fc Lissauer . 200d : 



Quintana et al. I . l2007r i. Figure 9.17 shows the results of some of their simulations. 



As shown here, depending on the mass-ratio of the binary and the initial values of 
its orbital parameters, in a few hundred million years, terrestrial planets can form 
around a close (0.01 to 0.1 AU) binary star system. 

Quintana and collea gues also studied terre strial planet formation in binaries 
with larger separations ( Quintana et al. I . 2007 ). Figure 9.18 shows the results of 
their simulations for a binary with a separation of 20 AU. Similar to figure 9.17, 
terrestrial-type objects are formed around the primary of the binary in a few hun- 
dred million years. Statistical analysis of their results, as shown in figures 9.19, 

9.20, and 9.21 indicate that, as expected in binaries with larger perihelia, where 
disk truncation has been smaller and more planet-forming material is available, 
terrestrial planet formation is elRcient and the numbe r of fin al terrestrial planets 
is large. The results of simulations by Quintana et al. ( 2007t ) also indicate that in 
a binary with a periastron distance larger than 10 AU, terrestrial planet formation 
can proceed efficiently in a region within 2 AU of the primary star. In binaries with 
periastra smaller than 5 AU, this region may be limited to inside 1 AU (figure 9.20). 
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Fig. 1.17. Terres t rial planet formation around a close binary system 
l|Quintana fc Lissauer I . |2006| ). The binary is circular and its separation is 0.05 AU. 
Each star of the binary has a mass of 0.5 solar-masses. A Jupiter-like planet has also been 
included in the simulation. The circles represent planetary embryos and planetesimals 
with radii that are proportional to their physical sizes. As shown here, the perturbative 
effect of the outer giant planet excite the orbits of the bodies at the outer edge of the 
disk and causes radial mixing as well as truncation. Within the first 100 Myr, several 
terrestrial-class objects are formed around the binary system. 



Despite the destructive role of the binary companion in increasing the relative 
velocities of planetesimals, which causes their collisions to result in erosion, this 
efficiency of terrestrial planet formation in binary systems may be attribute to the 
fact that the effect of the binary companion on increasing the relative velocities of 
planetesimals can be counterbala nced by dissipativc forces such as gas-drag and dy- 
namical friction ( Marzari et al. 1 . 1 199 7.: Marzari & SchoU , _200Q,) . The combination 
of the drag force of the gas and the gravitational force of the secondary star may 
result in the alignment of the periastra of planetesimals and increases the efficiency 
of their accretion by reducing their relative velocities. This is a process that is more 
effective when the sizes of the two colliding planetesimals are comparable and small. 
For colliding bodies with different sizes, depending on the size distribution of small 
objects, and the radius of each individual planetesimal, the process of the alignment 
of periastra may instead increase the relative velocities of the two objects, and cause 
their collisions to become eroding (figure 9.22, also see Thebault et al. . 2006h . 



22 Nader Haghighipour 







START 


1b = 


5 AU 


1b = 


7.5 AU 


q, = 10 AU 




0.3 
0.2 
















0.1 
0.0 
























5 


U.s: MYK 






O.fd MYK 




n Q \/VD 
U.(C MIK 


















U 


0.3 
0.2 
0.1 
0.0 


A; 






v: .■■■■>.v-.;-- . 
. y:«;.v ^, ■ 










10 MYR 






10 MYR 




1 MYR 


















tu 


0.3 
0.2 
0.1 
0.0 


• • 
















20 MYR 






20 MYR 




lOU MIK 


















QJ 


0.3 
0.2 
0.1 
0.0 


• 

• 






• 

• 

• 




• • ** 




0.5 


■ 100 MYR 






100 MYR 




100 MYR 




0.4 














tD 


0.3 
2 
0.1 
0.0 


• 

■ 






• 

• 




• 

• 




0.5 


200 MYR 






200 MYR 




200 MYR 




0.4 
















0.3 
0.2 












• 




0.1 
0.0 








O 

•° * 




: « o 




0.5 1.0 

a 


1.5 

(AU) 


2.0 


0.5 1.0 1.5 
a (AU) 


2.0 


0.5 1.0 1.5 2.0 
a (AU) 



Fig. 1.18. Terr e strial planet formation around the primary of a binary star system 
l|Quintana et al.ll2007^ . The stars of the binary are 0.5 solar- mass with semimajor axis 
of 20 AU. The eccentricity of the binary is 0.75 (left column), 0.625 (middle column), and 
0.5 (right column). As shown here, in each simulations, two terrestrial-type objects are 
formed after 100 Myr. The last row shows the results of additional simulations of the same 
systems, with final results showing in black, gray and white, corresponding to different 
runs. The differences in the final assembly of the planetary system of each simulation are 
results of the stochasticity of this type of numerical integrations. 



Simulations of terrestrial planet formation have also been extended to binaries 



with larger separations (20-40 AU) that also host a giant planet (jHaghighipour fc Raymond 



l2007l) . As discussed in the next section, by numerically integrating the orbits of the 
binary, its giant planet, and a few hundred planetary embryos, these authors have 
shown that it is possible to form Earth-like objects, with considerable amount of 
water, in the habitable zone of the primary of a moderately close binary-planetary 
system. 
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Fig. 1.19. Graphs of the final masses of the terrestrial planets formed in systems of figure 
9.18. The simulations have been run for three different masses of the binary stars. The 
red corresponds to simulations in a binary with 0.5 solar-masses stars, the blue represents 
results in a binary with 1 solar-mass stars, yellow is for a binary with a 1 solar-mass 
primary and a 0.5 solar-masses secondary, and black represents a binary with a 0.5 solar- 
masses primary and a 1 solar-mass secondary. These results show that despite the disk 
truncation in binaries wi th smaller perihe l ia, th e average masses of the final planets are 
not significantly altered. (jQuintana et al.ll2007l ). 
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Fig. 1.20. Graphs of the semimajor axis of the outermost planet of the simulations of 
figure 9.18. A shown here, while the outer edge of the disk is affected by the presence of 
the binary companion, the inne r portion of the disk, w here terrestrial planets are formed, 
stays unaffected by this object (jQuintana et al. U2007I '). 
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Fig. 1.21. Number of final terrestrial-type planets formed in the binaries of figure 9.18. 
As expected, for a given binary mass-ratio, the number of terrestrial planets increases in 
systems with larger p erihelia. This number a lso increases when the mass of the binary 
companion is smaller. l|Quintana et al.] ■ l2007^ . 
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Fig. 1.22. Encounter velocities of planetesimals with different sizes in a binary system wit h 
semimajor axis of 10 AU, eccentricity of 0.3, and mass-ratio of 0.5 l|Thebault et al. Il2006l ). 
The simulations include gas-drag. The vertical line shows the time of orbital crossing. As 
shown here, gas-drag lowers the encounter velocities of smaller equal-size planetesimals 
through the periastron alignment process. 
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1.4 Habitability 

It is widely accepted that a planet capable of developing life (as we know it), has to 
be able to continuously maintain liquid water on its surface and in its atmosphere. 
The capability of a planet in retaining water depends on its size and the processes 
involving its interior dynamics and atmospheric circulation. It also depends on its 
orbital parameters (i.e., its semimajor axis and orbital eccentricity) and the bright- 
ness of the central star at the location of the planet. While a dynamic interior and 
atmospheric circulation are necessary for a habitable planet to develop CO2 cycle 
and generate greenhouse effect (which helps the planet to maintain a uniform tem- 
perature), a long-term stable orbit, at a right distance from the star, is essential 
to ensure that the planet will receive the amount of radiation that enables it to 
maintain liquid water on its surface. These seemingly unrelated characteristics of 
a potential habitable planet, have strong intrinsic correlations, and combined with 
the luminosity of the star, determine the system's habitable zone. 

The inner and outer boundaries of a habitable zone va ry with the star's lu- 



minos i ty and the planet's atmospher i c circ ulati on models [seelMenou fc Tabachnik 



(|2003f) . |Jones. Underwood, fc Sleep I (|2005[ ). and lJones. Sleep, fc Underwood I (|2006( ) 



for a table of distances of the inner and outer boundaries of the habitable zones of 
exoplanetary systems]. A conservative estimate of the habitable zone of a star can 
be made by assuming that its inner edge is located at a distance closer than which 

water on the surface of the planet evaporates due to runaway greenhouse effect, 
and its outer edge is at a distance where, in the absence of CO2 clouds, runaway 
glaciation will freez e the water and creates permanent ice on the surface of the 
planet. As shown bv lKasting. Whitmire. fc Reynolds" (1993), such a definition of a 



habitable zone results in a habitable region between 0.95 AU and 1.15 AU for the Sun 
(figure 9.23). Th is is a somewhat conservative estim ate of the Sun's habitable zone 
and as noted bv I Jones. Un derwood, fc Sleep' (2005*), the outer edge of this re gion 
may, in fact, be farther away (Forget fc Pierrehumbert , , 1997.; . Williaras fc Kasting I . 
1997llMischna et al. Llioool i. 



Since the notion of habitability is based on life on Earth, it is possible to calculate 
the location of the boundaries of the habitable zone of a star by comparing its 
luminosity with that of the Sun. For a star with the surface temperature T and 
radius R, the luminosity L is given by 

L{R,T) ^ AnaT^R^, (1.5) 

where a is Boltzmann constant. Using equation (9.5) and the fact that Earth is in 
the habitable zone of the Sun, the radial distances of the inner and outer edges of 
the habitable zone of a star can be obtained from (Haghighipour 2006) 

The quantities Tq and i?© in equation (9.6) are the surface temperature and radius 
of the Sun, respectively, and r© represent the radial distance of Earth from the 
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Fig. 1.23. Habitable zone (|Kasting. Whitmire. fc Reynolds], 119931 ). 



Sun. Equation (9.6) implies that a habitable zone can be defined as a region around 
a star where an Earth-like planet can receive the same amount of radiation as 
Earth receives from the Sun, so that it can develop and maintain similar habitable 
conditions as those on Earth. 

As mentioned above, the orbit of a potential habitable planet in the habitable 
zone of a star has to be stable over long durations of time. As shown in section 9.2, 
the stability of the orbit of a planet in a binary system is strongly affected by the 
orbital motion of the binary companion. In binary systems where the primary hosts 
other planetary bodies (e.g., giant planets), the dynamics of a habitable planet will 
also be affected by the gravitational perturbations of these objects. It is therefore 
important to determine under what conditions a terrestrial-class object will have a 
long-term stable orbit in the habitable zone of a binary system, prior to construction 
a theory for the formation of Earth-like planets in such environments. 

Since a terrestrial-class planet is approximately two orders of magnitude less 
massive than a Jovian- type object, it will not have significant effect on the motion 
of the stars and the giant planets of a binary system. Therefore, as explained in 
sections 9.2.1 and 9.2.2, if a binary system does not contain a Jupiter-like planet, 
any dynamical criterion that is obtained for the stability or instability of a general 
planetary body, can also be applied to the dynamics of a terrestrial planet. Equation 
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Fig. 1.24. Graph of the hfetime of an Earth-size object in a circular orbit around the 
primary of 7 Cephei. The habitable zone of the primary has been indicated by HZ. No 
planet was placed in the region between the aphelion and perihelion distances of the giant 
planet of the system. As shown her e, only Earth-size pl anets close to the primary star 
maintain their orbits for long times (jHaghighipour I . I2OO6I ). 



(9.1) and the stability conditions presented by figure 9.5 can be used to determine 
the long-term stability of an Earth-like object in a binary system. 

If a binary contains giant planets, however, the situation is different. The grav- 
itational perturbations of the latter objects will have significant effects on the mo- 
tion and dynamics of terrestrial planets in the system. As shown by Haghighipour 
(2006), an Earth-size object, in a region between the giant planet and the primary of 
7 Cephei binary system, can maintain its stability only in orbits close to the primary 
star and outside the influence zoneH of the giant body. Integrating the equations of 
motion of a full four-body system, this author has shown that an Earth-like planet 
will not be able to sustain a stable orbit in the habitable zone of 7 Cephei's primary 
star (figure 9.24). However, it is possible for such an object to have a stable orbit 
when 0.3 < < 0.8 AU, 0° < ir = ip < 10°, and et < 0.4. Here qt represents 
the semimajor axis of the terrestrial planet and is its orbital inclinations with 
respect to the plane of the binary. 

As mentioned above, the instability of an Earth-like planet in the habitable zone 
of 7 Cephei can be attributed to the interaction between this object and the giant 
planet of the system. When the Earth- like planet is outside the giant planet's in- 

* The influence zone of a planetary object with a mass around a star with a mass M 
is deflned as the region between 3Rh — ap{l — Cp) and 3Rh + ap{l + ep) , where ap is 
the semimajor axis of the planet, and Rh ~ ap{mp/ZMY^^ is its Hill radius 
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Fig. 1.25. Formation of Earth-like planets in a binary-planetary systems. The top panel 
shows simulations in a binary with a separation of 30 AU, eccentricity of 0.2, and stellar 
components of 1 solar-mass. As shown here, an Earth-like planet (1.17 Earth-masses) with 
a water to mass ratio of 0.00164, is formed in the habitable zone of the primary at 1.16 AU, 
with an eccentricity of 0.02. The bottom panel shows the formation of an Earth-like object 
in a binary with a solar-mass primary and a 1.5 solar-masses secondary. The separation of 
the binary in this case is 30 AU, the mass of the Earth-like planets is 0.95 Earth-masses 
and its water to mass ratio is 0.00226. The semimajor axis of this planet and its orbital 
eccentricity are equal to 0.99 AU and 0.07, respectively. For the sake of comparison. Sun's 
habitable zone is approximately at 0.95-1.15 A U, Earth's orbital eccentricity is 0.017, and 
Earth's water to mass ratio is ~ 0.001 (|Haghighipour fc Raymond , . .2007. 1 . 
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Fig. 1.26. Formation of Earth-like planets in diflerent binary-planetary systems. As 
shown here, for a given binary mass-ratio, the delivery of water to terrestrial re- 
gion becomes less efficient as th e perihelion distance of the binary becomes smaller 
ijifaghighipour fc Ravmond I boO?! ') . 



fluence zone (e.g., at closer distances to the primray star) it can maintain its orbit 
for several hundred million years. Figure 9.24 suggest that, in order for a binary- 
planetary system to be habitable, its habitable zone has to be outside the influence 
region of its giant planet. In an S-type binary-planetar y system, this implies a pri- 



mary with a close-in habitable region. In a recent article. lHaghighipour fc Raymond 



(|2007( ) have studies the habitability of such a system. By considering a binary with a 
Sun-like primary star and a Jupiter-sized plane t in a circular orbit at 5 AU, and by 
adopting the model of iMorbidelli et al~ ( 2000( ). which is based on the assumption 
that water-carrying objects, in the Sun's asteroid belt, were the primary source of 
the delivery of water to Earth, these authors integrated the orbits of a few hundred 
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Fig. 1.27. The (e;,, a-^) parameter-space of an equal-mass binary-planetary system. Circles 
correspond to binaries with initial parameters chosen from figure 9.25, in which habitable 
planets were formed. Triangles represent systems in which the giant planet became unsta- 
ble. The numbers associated with each circle represents the mean eccentricity of the giant 
planet of the system at the end of the simulation. As shown here, moderately close bina- 
ries with lower eccen tricities (larger perihelia) are mor e suitable places for the formation 
of habitable planets l|Haghighipour fc Raymond 1 . 1200^ . 



protoplanetary (Moon- to Mars-sized) objects, and showed that it is indeed possible 
to form Earth-sized planets, with substantial amount of water, in the habitable zone 
of the primary star (figure 9.25). As shown by these authors, the mass and orbital 
parameters of the secondary star play important roles in the radial mixing of pro- 
toplanetary objects and the delivery of water to the habitable zone of the primary 
star. The giant planet of the system also plays the important role of transferring an- 
gular momentum from the secondary star to the disk of protoplanets, and enhancing 
the radial mixing of these objects. As shown in figure 9.26, water delivery is less 
efficient in binaries with smaller perihelia since in such systems, the close approach 
of the binary companion to the giant planet increases its eccentricity, which in turn 
results in stronger interaction between this object and the disk of protoplanets, 
ca using them to become unstable in v ery short time. The results of the simulations 
by iHaghighipour fc Raymond ( 2007 ) indicate that binary-planetary systems with 
giant planets at 5-10 AU, and binary perihelion distances of approximately 20 AU 
to 25 AU, will be more efficient in forming and hosting habitable planets (figure 
9.27). 
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1.5 Future Prospects 

The discovery of planets in binary star systems is one the interesting surprises of 
modern astronomy. Despite the long history of the study of planets in such environ- 
ments, the recent detection of planets in moderately close binaries has confronted 
astronomers with many new challenges. Many aspects of the formation process of 
these planets are still unresolved, and questions regarding their frequency and de- 
tection techniques demand more detailed investigation. 

The habitability of binary systems is also an open question. Although recent 
simulations of the late stage of terrestrial planet formation in binary systems have 
indicated that water-carrying planets can form in the habitable zone of a binary- 
planetary system, more studies are necessary to understand how protoplanetary 
objects can develop and evolve in such an environment. Such studies have impli- 
cations for investigating the habitability of extrasolar planets, and tic directly to 
several of near-future NASA missions, in particular, the space mission Kepler. 
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Fig. 1.28. Artistic rendition of tlie view from the moon of a giant planet in a triple star 
system. Figure from JPL-Caltech/NASA. 



